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This  paper  presents  a  dynamic  model  for  studying  the  cold  start  dynamics  and  observer  design  of  an 
auxiliary  power  unit  (APU)  for  automotive  applications.  The  APU  is  embedded  with  a  solid  oxide  fuel  cell 
(SOFC)  stack  which  is  a  quiet  and  pollutant-free  electric  generator;  however,  it  suffers  from  slow  start 
problem  from  ambient  conditions.  The  SOFC  APU  system  equips  with  an  after-burner  to  accelerate  the 
start-up  transient  in  this  research.  The  combustion  chamber  burns  the  residual  fuel  (and  air)  left  from 
the  SOFC  to  raise  the  exhaust  temperature  to  preheat  the  SOFC  stack  through  an  energy  recovery  unit. 
Since  thermal  effect  is  the  dominant  factor  that  influences  the  SOFC  transient  and  steady  performance,  a 
nonlinear  real-time  sliding  observer  for  stack  temperature  was  implemented  into  the  system  dynamics  to 
monitor  the  temperature  variation  for  future  controller  design.  The  simulation  results  show  that  a  100  W 
APU  system  in  this  research  takes  about  2  min  (in  theory)  for  start-up  without  considering  the  thermal 
limitation  of  the  cell  fracture. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Since  more  and  more  electronic  devices  are  integrated  into  the 
vehicles  as  the  standard  equipment,  the  demand  of  electricity  has 
approached  to  the  limit  of  the  current  energy  storage  batteries. 
An  independent  power  generation  system  with  compact  size  and 
excellent  efficiency  is  in  need.  Fuel  cells  are  one  of  the  solutions 
to  current  automotive  designers.  Delphi  and  BWM  have  developed 
a  solid  oxide  fuel  cell  (SOFC)  as  an  auxiliary  power  unit  (APU)  for 
transportation  applications  since  2000  [1-3].  The  SOFC  is  a  high 
temperature  operation  device  with  excellent  thermal  efficiency 
may  exceed  50%.  It  can  be  fueled  by  the  same  diesel  or  gasoline 
used  by  the  internal  combustion  engine,  or  even  multi-fuel  option, 
provided  that  an  external  reformer  is  installed.  The  tolerance  of 
hydrogen  impurity  also  makes  the  SOFC  a  suitable  choice  to  be 
the  automotive  APU.  Flowever,  the  start-up  for  this  kind  of  high 
temperature  fuel  cell  from  ambient  is  still  a  barrier  for  the  automo¬ 
tive  application.  In  this  paper,  the  high  temperature  exhaust  from 
a  burner  in  the  APU  system  was  used  to  accelerate  the  heating  up 
of  the  fuel,  the  air  and  the  SOFC  stack  to  reduce  the  cold  start  lag. 
A  temperature  sliding  observer  was  designed  to  monitor  the  tem¬ 
perature  variation  of  every  important  component  based  on  some 
detectable  system  parameters. 
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2.  SOFC  APU  system 

Fig.  1  shows  the  schematic  hardware  design  of  the  SOFC  APU 
system  in  this  research.  Inputs  of  the  SOFC  are  the  compressed 
air  and  fuel  (97%  H2  and  3%  water  vapor).  The  energy  recovery 
unit  (ERU)  includes  two  plate  heat  exchangers  to  recover  the  ther¬ 
mal  energy  from  the  high  temperature  exhaust  to  preheat  the 
inlet  air  and  the  inlet  fuel  of  the  SOFC  stack.  A  controller  was 
designed  to  regulate  the  valve  flows  via  valve  #1  and  #2  in  the 
diagram.  In  this  paper,  a  dynamic  model  of  the  SOFC  APU  sys¬ 
tem  was  developed  based  on  the  mass  and  energy  conservation 
laws  in  the  system.  The  dynamic  system  model  was  coded  on  a 
Matlab/Simulink  platform  for  on-line  real-time  simulation.  A  non¬ 
linear  sliding  observer  was  then  developed  and  embedded  into 
the  software  SOFC  system  to  carry  out  the  temperature  moni¬ 
toring  instantaneously.  In  this  paper,  the  SOFC  stack,  the  after 
-burner,  the  ERU  and  the  nonlinear  observer  are  the  major  sub¬ 
systems.  Their  dynamics  are  described  in  the  following  sections. 

2.1.  Solid  oxide  fuel  cell 

Solid  oxide  fuel  cells  are  one  of  the  high  temperature  power 
generators.  The  high  operating  temperature  enables  the  SOFC  to 
achieve  a  high  thermal  efficiency  with  excellent  tolerance  of  fuel 
impurity,  comparing  with  other  types  of  low  temperature  fuel  cells. 
The  SOFC  model  in  this  paper  includes  two  dynamic  systems:  the 
electrochemical  system  and  the  thermo-fluid  system.  The  former  is 
a  fast  response  reaction  system,  provided  that  the  temperature  is 


518 


P.-H.  Lin,  C.-W.  Hong /Journal  of  Power  Sources  187  (2009)  517-526 


Nomenclature 

A  area(m2) 

cw  heat  capacity  of  the  wall  (kj  kg-1 1<-1 ) 

E  output  voltage  of  solid  oxide  fuel  cell  (V) 

Eideai  ideal  voltage  of  solid  oxide  fuel  cell  (V) 

F  Faraday  constant  (C  mol-1 ) 

Ft  correction  factor  of  the  heat  exchanger 

/  range  of  parameter  uncertainty 

i  current  density  (A  cm-2 ) 

i0  exchange  current  density  (A  cm-2 ) 

inmit  limiting  current  density  (A  cm-2 ) 

I  current  demand  (A) 

h  enthalpy  (kj  kg-1) 

ho  stagnation  enthalpy  (kj  kg-1 ) 

Nc  number  of  the  channels  of  the  heat  exchanger 
n  number  of  transfer  electrons 

m  mass  (kg) 

/<i ,  k2  sliding  gain 

P  pressure  (kPa) 

P  estimated  pressure  (kPa) 

Q  heat  (J) 

R  gas  constant  (J  mol-1  K-1 ) 

Rohm  fuel  cell  ohmic  resistance  (£2) 

S  entropy  (JK_1) 

T  temperature  (V) 

T  estimated  temperature  (K) 

A Tim  logarithmic  mean  temperature  different 

t  elapsed  time  (s) 

u  specific  internal  energy  (kj  kg-1 ) 

U  overall  heat  transfer  coefficient 

Greek  letters 

a  charge  transfer  coefficient 

rjact  activation  overpotential  (V) 

Tjcon  concentration  overpotential  (V) 

pohm  ohmic  overpotential  (V) 

Subscripts 

HE  heat  exchanger 

amb  ambient 

an  anode 

ca  cathode 

cold  cold 

comb  combustion 

8 

hot  hot 

in  inlet 

j  species  of  inlet  gas 

PEN  positive-electrolyte-negative 

out  outlet 

w  wall 


raised  to  the  operation  criterion  (800 1<  approximately).  The  latter 
is  much  slower  because  it  involves  the  flow  inertia  and  thermal 
inertia  in  the  transport  system,  which  consists  of  the  pipes,  valves, 
damping  volumes,  and  heat  exchangers  in  Fig.  1.  The  time  lag  of 
fluid  dynamics  and  heat  transfer  among  the  flow  transport  is  the 
major  concern  in  the  design  of  this  fast  start-up  SOFC  system. 

The  thermo-fluid  system  is  modeled  by  a  modified  filling-and- 
emptying  technique  which  was  developed  for  the  turbocharged 
engine  simulation  originally  [4].  Both  temperature  and  pressure 
dynamics  are  derived  from  the  mass  and  the  energy  conservation 
laws  plus  the  differential  form  of  the  gas  state  equation.  The  energy 


(temperature)  and  mass  dynamics  can  be  expressed  by 


J^jdQldt)  +  hoj(dmjldt )  -  u(dm/dt )  ]  (1/m) 


dT 

dt 


du/dT 


(1) 

(2) 


where  ^^( dmj/dt )  is  the  summation  of  the  mass  that  enters  and 

j 

leaves  the  control  volume.  The  term  ^^(dQ/dt)  includes  all  the 

heat  transfer  (through  the  control  surface)  and  additional  heat 
release  (due  to  combustion  or  electrochemical  reaction  inside  the 
control  volume);  h0j{dmj/dt )  is  the  summation  of  the  change  of 

j 

the  specific  stagnation  enthalpy  due  to  the  mass  flow  in  and  flow 
out  of  the  control  volume;  u  dm\dt  is  the  internal  energy  variation 
due  to  the  system  mass  change;  du/dT  is  the  specific  heat  at  con¬ 
stant  volume.  The  pressure  dynamics  inside  the  control  volume  can 
then  be  derived  from  the  state  equation  of  ideal  gas,  which  takes 
the  form  of 


dP  _  (  \  dm  1  dT\ 

dt  ~  \m  dt  +  T  dt ) 


(3) 


The  electrochemical  reaction  rate  is  much  faster  than  the  flow  of 
the  thermo-fluid  system.  Hence,  the  output  voltage  of  the  electro¬ 
chemical  system  of  the  SOFC  can  be  evaluated  quasi-statically  and 
is  expressed  by 


E  —  Eideal  ~  Pact  ~  Peon  ~  Pohm  (4) 

where  Eideai  is  the  theoretical  Nernst  voltage  of  the  SOFC;  pact , 
Peon  and  pohm  are  the  activation  overpotential,  the  concentra¬ 
tion  overpotential  and  the  ohmic  overpotential,  respectively.  These 
overpotentials  are  related  with  the  operation  current  density,  z,  and 
the  instantaneous  stack  temperature,  T,  as  below 

RT  RT 

Pact  =  a  b\n  i  where  a  « - —  In  z„  and  b^ — -  (5) 

naF  naF 


Pohm  ~  I  Rohm  (7) 

where  a  and  b  are  the  span  and  the  slope  of  the  Tafel  plot  of  the 
activation  overpotential;  n  is  the  number  of  transfer  electrons;  a  is 
the  charge  transfer  coefficient;  F  is  the  Faraday  constant;  R  is  the 
gas  constant;  T  is  the  operation  temperature;  z'o  is  the  exchange 
current  density;  is  the  limiting  current  density;  Rohm  is  the 
total  resistance  of  the  electrolyte,  and  /  is  the  current  demand.  All 
the  terms  in  the  right-hand-side  of  the  above  equations  are  func¬ 
tions  of  the  stack  temperature  T,  and  temperature  is  a  function  of 
time  in  our  case.  Hence  the  first-order  nonlinear  differential  equa¬ 
tions  in  Eqs.  (1-3)  configure  the  system  dynamics  of  the  control 
volume.  The  SOFC  stack  temperatures,  including  the  anode  cham¬ 
ber  temperature,  cathode  chamber  temperature,  and  the  electrolyte 
(or  positive-electrolyte-negative,  PEN)  temperature,  are  the  major 
parameters  to  be  observed  for  control  purposes  [5,6]. 

An  orifice  analogy  model  to  simulate  the  mass  flow  across  the 
boundary  due  to  pressure  difference  between  control  volumes  has 
been  employed.  The  mass  flow  rate  is  calculated  by  an  isentropic 
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Fig.  1.  Schematic  of  a  typical  SOFC  APU  system. 


Fig.  2.  Block  diagram  of  the  SOFC  stack  dynamics  calculation. 


nozzle  flow  multiplied  by  a  discharge  coefficient  as 


(8) 


where  Q  is  the  discharge  coefficient;  Ab  is  the  flow  area  at  the 
boundary;  Pu,  Pd,  Pr  are  the  upstream,  downstream  pressures  and 
the  pressure  ratio;  Tu  is  the  upstream  temperature;  R  is  the  gas 
constant  and  y  is  the  specific  heat  ratio. 

Fig.  2  shows  the  block  diagram  of  the  calculation  of  the  SOFC 
stack  dynamics.  The  SOFC  model  includes  the  blocks  of  the  anode 
dynamics,  the  cathode  dynamics  and  the  PEN  dynamics.  The  upper 
left  side  of  the  figure  shows  the  detailed  flow  diagram  of  the  calcu¬ 
lation  of  the  anode  dynamics.  The  cathode  dynamics  resemble  the 
anode  dynamics.  Note  that  the  inputs  of  the  whole  block  diagram 
are;  Pamb .  Tamb,  Ph2  .  Th2  .  Pair ,  Pair  and  the  current  demand  /;  the  out- 
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Table  1 

The  specification  and  operating  conditions  of  the  SOFC  APU  system. 


Component 

Specification  and  operating  conditions 

Solid  oxide  fuel  cell 

Type:  Co-flow  planar  SOFC 

Stack  volume:  10  cm  x  25  cm  x  15  cm 
Single  cell  reaction  area:  50  cm2 

Numbers  of  cells:  5 

Suggested  operation  temperature:  1100K 
Fuel:  97%  hydrogen  +  3%  water  vapor 

Input  pressure:  150  kPa 

Input  temperature:  300  K 

Oxidizer:  air 

Input  pressure:  180  kPa 

Input  temperature:  300  K 

Rated  power:  100  W 

Heat  exchanger 

Plate  type 

Channel  numbers:  16 

Heat  transfer  coefficient:  929  Wm~2 1<-1 
Exchange  area:  0.0729  m2 

Correction  factor:  1 

put  variables  ate  the  stack  voltage  E  and  the  PEN  temperature  TPEN. 
A  3%  of  water  vapor  is  added  into  the  H2  inlet  stream  for  protec¬ 
tion  purpose.  The  H2  consumption  rate  can  be  estimated  from  the 
Faraday  law,  that  is 

IMh2 

"%=  nF 

O) 

where  I  is  the  current  demand  from  the  external  load;  M  is  the 
molecular  weight;  n  is  the  number  of  transfer  electrons;  F  is  the 
Faraday  constant. 


where  m  is  the  consumption  (or  production)  rate,  h  is  the  enthalpy. 
Subscripts/,  p,  r  denote  formation,  product  (water  vapor)  and  reac¬ 
tants  (hydrogen  and  oxygen);  superscript  0  means  at  standard 
condition  (1  atm,  298 1<).  Note  that  an  iteration  process  is  required 
to  predict  the  burnt  (product)  temperature  and  the  heat  release  in 
the  same  time  (as  described  by  typical  thermodynamics  textbooks). 

2.3.  Energy  recovery  unit 

The  energy  recovery  unit  (ERU)  retrieves  heat  from  the  high  tem¬ 
perature  exhaust  of  the  burner  and  transfers  it  to  pre-heat  the  SOFC 
stack.  It  consists  of  two  plate  heat  exchangers  which  are  used  to 
heat  up  the  incoming  fuel  and  the  incoming  air  respectively.  The 
heat  transfer  rate  of  the  plate  heat  exchanger  can  be  expressed  by 
[7] 

=  (Nc  -\)AHEUFTATlm  (11) 

where  Qhe  is  the  heat  transfer  rate;  Nc  is  the  number  of  channels; 
Ahe  is  the  effective  heat  exchange  area;  U  is  the  overall  heat  transfer 
coefficient;  FT  is  its  correction  factor  [8] ;  and  A T/m  is  the  logarithmic 
mean  temperature  difference,  which  is  defined  by 

_  ( Thot,in  ~  Tcold.out )  —  i^hot,out  ~  ? cold, in )  (12) 

hot, in  ~  Tcold.out)/ (Thot,out  ~  T cold, in )] 

The  above  equations  are  used  to  evaluate  the  heat  transfer  rate  and 
the  outlet  gas  temperature  at  steady  state  operation. 

2.4.  Transient  wall  temperature  evaluation 


2.2.  After-burner 

An  after-burner  is  placed  at  the  exit  of  the  SOFC  stack  to  raise 
the  exhaust  gas  temperature  by  burning  the  residual  fuel.  This  high 
temperature  exhaust  gas  is  then  fed  back  as  the  heat  source  for 
those  two  heat  exchangers  in  the  energy  recovery  unit.  The  burner 
is  considered  as  an  adiabatic  system  so  that  the  heat  released  from 
combustion  is  totally  utilized  to  increase  the  gas  temperature.  The 
temperature  and  the  pressure  dynamics  of  the  burner  are  calcu¬ 
lated  from  Eqs.  (1 )  and  (3).  An  additional  term  dQcombldt  was  added 
in  Eq.  (1)  to  take  account  of  the  chemical  energy  release  during 
combustion. 

The  combustion  heat  release  rate  is  derived  from  the  enthalpy 
change  rate  between  the  product  (at  product  temperature  Tp)  and 
the  reactants  (at  reactant  temperature  Tr).  It  is  expressed  by 

=  mp[Ah°fp  +  hp(Tp)  -  h°p]  -  JfmAA  h°fr  +  hr(Tr )  -  h°r] 

(10) 


The  transient  wall  temperature  (such  as  the  PEN  temperature) 
is  evaluated  from  a  first-order  heat  transfer  equation  as  below: 

mwcw  dt  =  hgAw(Tg  —  Tw)  (13) 

where  mw  is  the  wall  mass,  cw  its  heat  capacity  of  the  wall,  h  is 
the  effective  heat  convection  coefficient,  Aw  is  the  wall  area,  Tg  is 
the  heating  gas  temperature  and  Tw  is  the  wall  temperature.  Under 
transient  conditions,  a  change  in  the  effective  gas  temperature  does 
not  result  in  an  instantaneous  change  of  the  stack  wall  temperature 
due  to  thermal  inertia  of  the  system.  The  actual  wall  temperature 
at  time  step  n  + 1  is 

Tw,n+ 1  —  lw,n  +  {Tw,n+ l,r=0  —  7w,n)(l  —  e  At/T)  (14) 

The  thermal  time  constant  r  is  calculated  from  an  electric  circuit 
analogy  by 

^  _  xwpwcw  (15) 

^  +  xw/kw+l/h2 


Table  2 

Numerical  values  for  the  simulation  parameters. 


Component 

Parameter 

Value 

Volume  of  the  anode  chamber  (Van) 

0.0002  m3 

Volume  of  the  cathode  chamber  (Van) 

0.0002  m3 

Mass  of  the  electrolyte 

0.0003  kg 

SOFC 

Heat  capacity  of  the  electrolyte 

1000 J  kg-1  K"1 

Charge  transfer  coefficient  (a) 

0.3 

Exchange  current  density  (i'o) 

0.1  A  cm-2 

Limiting  current  density 

0.9  A  cm-2 

Ohmic  resistance  ( Rohm ) 

0.01  a 

After-burner 

Volume  of  the  burner  ( Vfaurner) 

0.002  m3 

Heat  exchanger 

Volume  of  the  heat  exchanger  chamber  (VHe) 

0.002  m3 

Sliding  gains  (k\,  k2) 

200,  20 

Observer 

Initial  estimated  pressure  (P) 

111.46  kPa 

Initial  estimated  temperature  (T) 

400  K 
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(a)  Current  demand  and  voltage  output 


Elapsed  Time  (s) 

Fig.  4.  Single  cell  voltage  output,  current  demand  and  PEN  temperature  versus  elapsed  time  showing  the  start-up  and  dynamic  performance  of  the  SOFC.  (a)  Current  demand 
and  voltage  output,  (b)  PEN  temperature. 


in  which,  xw  is  the  wall  thickness,  pw  is  the  wall  density,  cw  is  the 
specific  heat  capacity  of  the  wall,  kw  is  the  conduction  coefficient 
of  the  wall,  h\  and  /12  are  the  heat  convection  coefficients  of  both 
sides  of  the  wall. 


2.5.  Sliding  observer 

The  SOFC  stack  is  a  high  temperature  operation  device  which 
needs  to  be  thermally  insulated  in  a  chamber.  The  thermal  effect 
plays  the  dominant  factor  that  influences  the  SOFC  performance. 
A  novel  sliding  observer  was  designed  to  obtain  the  immeasur¬ 
able  gas  temperatures  inside  the  SOFC  system  via  some  measurable 
dynamic  signals  (such  as  pressure  signals  from  the  simulation). 
This  is  for  future  start-up  control  in  the  next  research.  Due  to  the 
fact  that  the  system  parameters  are  time  dependent,  the  observer 
must  be  also  a  temporal  function  and  possesses  the  robustness  for 
parameter  uncertainties.  The  sliding  observer  is  designed  based  on 
a  sliding  surface  concept,  in  which  a  switching  law  is  added  into 


the  state  dynamic  equations  of  the  whole  system.  The  switching 
law  is  a  function  of  the  error  between  the  target  system  state  and 
the  estimated  system  state.  This  observer  shows  sliding  behavior 
in  the  estimated  state  trajectory  as  shown  in  Fig.  3.  The  sliding 
behavior  enables  the  observer  estimate  the  actual  system  state  suc¬ 
cessfully  in  a  short  period  even  with  a  large  model  uncertainty  and 
measurement  disturbance  [9-12]. 

In  our  case,  we  use  the  observer  to  acquire  the  system  states 
that  are  difficult  to  be  obtained  in  the  realistic  condition,  such  as  the 
SOFC  gas  temperatures  in  the  example  system.  The  sliding  observer 
for  the  anode  dynamics  is  derived  below  as  an  example  for  the 
other  parts  of  the  SOFC  APU  system.  The  observer  is  built  up  to 
acquire  the  anode  gas  temperature  based  on  available  pressures 
whose  equations  are  expressed  by 


Tan  — 


Pan  =  Pan(drhan/(rhandt)  +  dfan/(fandt))  +  /<i  Sgn(Pan) 
dfhj 


dm;  ^dman  \  1 

.LdF+2>^r-u^r 


/^+k2sgn(Pan)  (16) 
d  Tan 


(U 

CL 


CD 

=5 
(/) 
0 0 
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(a)  Anode  side  heat  exchanger  pressure 


(b)  Anode  side  heat  exchanger  temperature 
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(c)  Cathode  side  heat  exchanger  pressure  (d)  Cathode  side  heat  exchanger  temperature 


Fig.  5.  Heat  exchanger  pressure  and  temperature  dynamics  of  the  fuel  (anode  side,  (a)  and  (b))  and  the  air  (cathode  side,  (c)  and  (d))  from  start-up. 
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(a)  Anode  pressure 


(b)  Anode  temperature 
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Fig.  6.  Dynamics  of  the  (a)  pressure,  (b)  temperature,  (c)  mass  and  (d)  energy  flow  of  the  anode  volume. 


where  pan  and  Tan  are  the  estimated  states;  k\  and  k2  are  the  slid¬ 
ing  gains  of  the  observer;  pan  is  the  error  between  the  estimated 
pressure  and  the  target  pressure.  The  sliding  states  are  set  to  be 
x  =  [  x\  x2  ]  =  [  Pan  T  ],  and  the  sliding  surface  is  chosen  as  the 
pressure  error  (xi )  which  is  the  difference  between  the  estimated 
pressure  and  the  target  pressure: 


The  state  trajectory  will  converge  to  s  =  0,  which  means  that  the 
error  x,\  will  be  zero,  if  the  following  equation  is  satisfied 


<  -fsgn(s) 


(19) 


S  =  X  1  -  X\  =X\ 

and 


( 17)  where  f  is  a  positive  constant  and  sgn  denotes  ‘sign’,  which  is  either 
positive  or  negative.  If  Axi  in  Eq.  (16)  is  bounded  by  a  constant/, 
then  the  range  of  the  sliding  gain  k\  is 


=  *i  -  xi  =  Axi  -  kisgn(xi ) 


(18)  /<!>/  +  ? 
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Fig.  7.  Dynamics  of  the  (a)  pressure,  (b)  temperature,  (c)  mass  and  (d)  energy  of  the  cathode  volume. 
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(a)  After-burner  pressure 


(b)  After-burner  temperature 
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Fig.  8.  Dynamics  of  the  (a)  pressure,  (b)  temperature,  (c)  mass  and  (d)  energy  of  the  after-burner. 


A  reduced  dimensional  sliding  observer  can  be  derived  for  x2 
from  Eq.  (15)  that 


~  d(Ax2)  _  k 2  d(A*i)  = 
2  dx2  /<!  dx2  - 


(21) 


where  A2,  B2 ,  C2  are  temporal  variables  and  A2  =  9(  Ax2)/9x2,  B2  = 
k2/k\  and  C  =  9(  Axi  )/9x2;  Ax2  is  the  difference  between  the  target 
state  x2  and  the  estimated  state  x2.  Eq.  (21 )  will  converge  to  zero  if 
a  proper  sliding  gain  k2  is  selected.  The  sliding  gain  k2  for  the  anode 
temperature  observer  can  be  obtained  by  the  following  equation: 


feiC2(jg-qg) 

A2 


(22) 


where  a  is  the  pole  of  A2,  and  is  the  new  pole  location.  Similar 
procedure  is  applied  to  the  cathode  observer,  and  the  representing 
PEN  temperature  is  then  calculated  based  on  the  estimated  anode 
and  cathode  temperatures. 


Gas  composition  in  the  after-burner 


Fig.  9.  Gas  composition  ratio  variation  inside  the  after-burner. 


3.  Simulation  results 

3.1.  Transient  performance  under  cold  start  condition 

The  dynamic  models  of  each  component  of  the  SOFC  APU  system 
were  built  up  and  coded  on  a  Matlab/Simulink  platform.  A  numer¬ 
ical  scheme  ODE4  (Runge-Kutta)  was  employed  to  solve  the  three 
dynamic  equations  for  mass,  temperature  and  pressure  respec¬ 
tively.  The  specification  of  the  system,  operating  conditions  and 
input  numerical  parameters  are  listed  in  Tables  1  and  2.  The  simu¬ 
lation  was  carried  out  on  a  typical  Pentium  4  personal  computer 
and  the  computation  time  is  approximately  half  of  the  realistic 
operation  time  of  the  SOFC  system. 

Fig.  4(a)  shows  the  current  demand,  single  cell  voltage  output 
and  the  electrolyte  temperature  variation  of  the  SOFC  stack  under 
a  cold-start  condition.  Two  stages  of  current  demand  were  given: 
10  A  current  was  applied  to  the  cell  when  the  electrolyte  tempera¬ 
ture  reached  800  K;  20  A  was  applied  at  the  300th  second  instantly 
to  test  the  cell  dynamics.  In  Fig.  4(a)  and  (b),  the  electrolyte  tem¬ 
perature  starts  from  300  K  and  reaches  800  K  at  the  23rd  second 
of  the  elapsed  time.  At  this  time  a  starting  current  10  A  is  applied. 
The  cell  voltage  in  Fig.  4(a)  increases  with  the  stack  temperature 
and  keeps  at  a  constant  output  at  about  0.8  V.  Fig.  4(b)  shows 
the  instantaneous  PEN  temperature  and  an  enlarged  view  of  the 
temperature  variation  is  displayed  below  to  trace  all  the  transient 
conditions.  The  PEN  temperature  starts  from  300 1<  and  keeps  at 
a  constant  temperature  (about  1095  K)  after  100  s.  The  magnified 
view  reveals  the  relation  between  the  temperature  variation  and 
the  current  command.  The  reason  for  the  temperature  variation 
is  due  to  the  heat  released  from  the  chemical  reaction.  The  heat 
release  becomes  greater  when  a  higher  current  demand  is  applied, 
so  that  the  temperature  climbs  to  a  higher  peak  proportionally  after 
300s. 

Fig.  5  shows  the  pressure  and  temperature  variation  of  the  pre¬ 
heated  hydrogen  stream  and  air  stream  inside  the  energy  recovery 
unit.  The  inlet  pressures  of  both  the  hydrogen  and  the  air  are 
assumed  to  be  150  kPa  and  180  kPa,  respectively;  the  inlet  tem¬ 
peratures  are  both  set  at  the  room  temperature.  The  compressed 
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Fig.  10.  Tracking  performance  of  the  sliding  observer:  (a)  anode  pressure;  (b)  cathode  pressure.  Enlarged  views  show  the  excellent  tracking  performance  of  the  observer. 


air  is  used  to  provide  more  oxygen  to  the  SOFC  and  the  burner. 
Fig.  5(a)  and  (c)  show  that  the  anode  and  cathode  pressures  vary 
at  the  beginning  and  then  keep  at  a  constant  value  near  142  kPa 
and  165  kPa  after  100  s.  They  are  both  lower  than  the  inlet  pres¬ 
sure  due  to  pressure  drops  across  boundaries.  Temperatures  of  the 
preheated  fuel  and  the  preheated  air  before  entering  the  fuel  cell 
increase  with  the  burner  temperature  and  reach  to  a  constant  value 
at  about  1028  K  and  1156  K,  respectively.  The  preheated  air  temper¬ 
ature  is  higher  than  the  fuel  because  of  the  difference  of  the  heat 
capacity  and  the  mass  in  the  system.  All  the  second  peaks  near 
the  300th  second  are  due  to  the  application  of  the  current  demand 
preset  at  20  A. 

The  dynamics  of  the  control  volumes  for  the  anode  and  the 
cathode  plenums  are  shown  in  Figs.  6  and  7,  respectively.  Pressure 
variation  in  the  beginning  of  the  simulation  can  be  observed  in  both 
anode  and  cathode  volumes.  This  is  due  to  the  imbalance  of  the 
mass  flow  rate  inside  the  channels  as  shown  in  Figs.  6(c)  and  7(c). 


Note  that  when  the  anode  pressure  reaches  steady  state,  the  net 
mass  flow  rate  of  hydrogen  and  water  vapor  in  the  anode  is  still 
imbalance  as  shown  in  Fig.  6(c).  This  is  due  to  an  additional  mass, 
oxygen  ions,  entered  the  anode  plenum  when  the  electrochemical 
reaction  started,  so  that  the  mass  flows  out  of  the  anode  becomes 
higher  in  order  to  stabilize  the  pressure.  In  Figs.  6(b)  and  7(b),  both 
temperatures  of  the  anode  and  the  cathode  rise  from  the  beginning 
and  stabilize  at  1032  K  and  1155  K,  respectively.  The  anode  tem¬ 
perature  and  the  cathode  temperature  are  higher  than  the  input 
temperature  because  an  additional  heat  release  from  the  electro¬ 


chemical  reaction,  i.e. 


in  the  energy  equation. 


The  values  of  the  net  energy  flow,  resulted  from  heat  transfer, 
enthalpy  flow  in  and  out,  internal  energy  change  and  the  electro¬ 
chemical  heat  release,  are  positive  in  Figs.  6(d)  and  7(d).  Hence,  it 
raises  the  gas  temperature  in  the  control  volumes. 
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Fig.  11.  Tracking  performance  of  the  sliding  observer:  (a)  and  (b)  anode  temperature;  (c)  and  (d)  cathode  temperature;  (e)  and  (f)  PEN  temperature.  Magnified  views  show 
the  estimated  temperature  (red  dashed  line)  catches  up  the  target  temperature  (solid  blue  line)  before  100  s  (within  2  min).  (For  interpretation  of  the  references  to  color  in 
this  figure  legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 


The  transient  performance  of  the  after-burner  is  shown  in  Fig.  8. 
The  sudden  increase  of  the  pressure  in  the  beginning,  shown  in 
Fig.  8(a),  is  due  to  the  combustion  heat  release  in  Fig.  8(d).  The 
burner  pressure  stabilizes  at  a  constant  value  near  104  kPa  after 
100  s.  The  increase  of  the  burner  temperature  in  Fig.  8(b)  is  also  due 
to  the  combustion  heat  release,  and  the  temperature  reaches  1160  K 
after  the  burner  is  ignited.  Fig.  8(c)  shows  that  the  mass  flow  in  and 
out  of  the  burner  balances  all  the  time,  except  in  the  beginning, 
which  means  that  the  mass  flow  in  an  open  system.  In  Fig.  8(d), 
combustion  heat  release  results  in  the  high  temperature  exhaust 
gas  at  the  burner  exit.  Fig.  9  shows  the  gas  composition  variation 
inside  the  burner.  During  the  unsteady  combustion  period,  the  ratio 
of  the  hydrogen  drops  from  0.88  to  0.79  and  then  returns  to  0.86. 
Due  to  the  drop  of  the  hydrogen  ratio,  the  other  gas  compositions 
slightly  increase  and  then  stabilize  at  a  higher  value  than  the  initial 
condition  after  100  s. 


3.2.  Tracking  performance  of  the  sliding  observer 

Fig.  10  shows  the  target  pressure  (solid  blue),  the  estimated  pres¬ 
sure  (dashed  red)  and  the  estimation  error  (dashed  green)  of  the 
anode  observer  and  the  cathode  observer.  In  Fig.  10(a),  the  anode 
observer  successfully  tracks  the  target  pressure  after  10  s.  The  mag¬ 
nified  view  for  the  error  shows  that  there  is  a  peak  which  is  about 
20  kPa  appearing  in  the  beginning  of  the  estimation  and  then  chat¬ 
ters  around  zero  during  all  the  other  transients  in  the  later  stage. 
Fig.  10(b)  shows  that  the  cathode  observer  successfully  tracks  the 
target  pressure  and  reduces  the  error  from  12  kPa  to  zero  in  70  s.  All 
the  other  transients  are  also  well-predicted. 

Fig.  11  shows  the  tracking  performance  of  the  temperature 
observer  at  the  anode,  the  cathode  and  the  PEN.  The  initial  tem¬ 
perature  for  the  observer  was  set  at  400  K  to  magnify  the  effect 
of  uncertainty,  while  the  actual  temperature  starts  from  300  K. 
Fig.  11(a),  (c),  and  (e)  indicate  that  the  sliding  observer  (dashed 
red)  tracks  the  actual  temperature  (solid  blue)  successfully  in  60  s. 
Fig.  11(b),  (d),  and  (f)  reveal  that  the  observed  anode  temperature, 
the  cathode  temperature,  and  the  PEN  temperature  take  approxi¬ 


mately  10,  75,  and  3  s  to  track  back  the  original  profile  due  to  the 
sudden  increase  of  the  current  at  the  300th  second.  All  the  above 
results  prove  that  the  sliding  observer  is  capable  of  feeding  back 
a  correct  temperature  to  the  SOFC  controller  which  is  still  under 
development. 

4.  Conclusions 

Cold  start  dynamics  of  a  100  W  SOFC  APU  system  have  been 
investigated  in  this  paper.  A  modified  filling-and-emptying  model 
was  employed  to  derive  the  dynamic  models  of  the  SOFC  sys¬ 
tem  based  on  the  mass  and  energy  conservation  laws.  The  high 
temperature  exhaust  of  the  embedded  burner  was  used  to  accel¬ 
erate  the  cold  start  of  the  APU  system  with  an  energy  recovery 
unit.  The  simulation  results  show  that  the  APU  system  tempera¬ 
ture  in  this  research  takes  100  s  in  theory  (without  considering  the 
material  thermal  stress  problem)  to  reach  the  steady  state.  A  non¬ 
linear  sliding  observer  was  designed  to  obtain  the  immeasurable 
gas  temperature  instantaneously  of  the  anode,  the  cathode  and  the 
electrolyte.  The  temperature  observer  is  able  to  acquire  the  accu¬ 
rate  state  value  in  a  short  period  with  a  proper  sliding  gain.  It  can 
be  concluded  at  this  stage  that  the  observer  is  useful  in  developing 
a  fast  start-up  controller  for  the  APU  system.  A  controller  using  the 
feedback  observer  temperature  to  track  the  gradient  of  the  assigned 
heating  up  profile  to  avoid  the  fracture  problem  is  under  develop¬ 
ment. 
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